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Transcriptomic characterization of two major Fusarium resistance 
quantitative trait loci (QTLs), Fhb1 and Qfhs.ifa-5A, identifies novel 
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SUMMARY 

Fusarium head blight, caused by Fusarium graminearum, is a 
devastating disease of wheat. We developed near-isogenic lines 
(NILs) differing in the two strongest known F. graminearum resist- 
ance quantitative trait loci (QTLs), Qfhs.ndsu-3BS (also known as 
resistance gene Fhb1) and Qfhs.ifaSA, which are located on the 
short arm of chromosome 3B and on chromosome 5A, respec- 
tively. These NILs showing different levels of resistance were used 
to identify transcripts that are changed significantly in a QTL- 
specific manner in response to the pathogen and between mock- 
inoculated samples. After inoculation with F. graminearum spores, 
16 transcripts showed a significantly different response for Fhb1 
and 352 for Qfhs.ifaSA. Notably, we identified a lipid transfer 
protein which is constitutively at least 50-fold more abundant in 
plants carrying the resistant allele of Qfhs.ifaSA. In addition to 
this candidate gene associated with Qfhs.ifaSA, we identified a 
uridine diphosphate (UDP)-glycosyltransferase gene, designated 
TaUGT12887, exhibiting a positive difference in response to the 
pathogen in lines harbouring both QTLs relative to lines carrying 
only the Qfhs.ifaSA resistance allele, suggesting Fhb1 depend- 
ence of this transcript. Yet, this dependence was observed only in 
the NIL with already higher basal resistance. The complete cDNA 
of TaUGT12887 was reconstituted from available wheat genomic 
sequences, and a synthetic recoded gene was expressed in a 
toxin-sensitive strain of Saccharomyces cerevisiae. This gene con- 
ferred deoxynivalenol resistance, albeit much weaker than that 
observed with the previously characterized barley HvUGT13248. 
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INTRODUCTION 

The plant-pathogenic fungus Fusarium graminearum frequently 
infects small grain cereals, particularly wheat (Leonard and 
Bushnell, 2003). Fusarium head blight (FHB) is one of the most 
important field crop diseases in temperate regions, resulting in 
severe yield losses and reduced grain quality by contamination 
with the trichothecene deoxynivalenol (DON; Desjardins, 2006). 
This mycotoxin is a serious threat to human and animal health 
(Pestka, 2010). Maximum tolerated levels in food have been 
enacted in the European Union (The European Commission, 2006), 
and advisory levels have been issued by the Food and Drug Admin- 
istration in the USA. As the control of FHB by crop rotation, tillage 
regimes or the use of fungicides is only partially successful, an 
understanding of the genetic control of FHB resistance and the 
deployment of resistance genes in wheat varieties are likely to be 
the most promising strategies. 

Genetic mapping of FHB resistance in wheat has led to the 
identification of several quantitative trait loci (QTLs) that confer 
partial resistance (Buerstmayr ef al., 2009). The two most effective 
QTLs exhibiting FHB resistance were identified in the Chinese 
Spring wheat cultivar Sumai3 and derivatives, and are positioned 
on the short arm of chromosome 3B {Fhbl) and on chromosome 
5A (QfhsJfaSA) (Anderson et al., 2001; Buerstmayr ef al., 2002, 
2003). Molecular markers have been developed for the accurate 
screening of Fhbl (Liu etai, 2008) and Qfhs.ifaSA (Buerstmayr 
ef al., 2002) resistance alleles. Fhbl and Qfhs.ifaSA employ dif- 
ferent mechanisms. Fhbl confers so-called type II resistance by 
slowing down or inhibiting the spread of the pathogen from the 
initial infection site, whereas Qfhs.ifaSA contributes to type I 
resistance by lowering the rate of initial infection and also confers 
type II resistance - albeit to a lesser extent. 

The genetic determinants underlying Fhbl and Qfhs.ifaSA are 
still unknown. Despite significant efforts, the Fhbl gene has not 
yet been isolated (Liu ef al., 2006, 2008). In addition, the release 
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of parts of the wheat chromosome 3BS sequence of the FHB- 
susceptible cultivar Chinese Spring, which should cover the Fhbi 
locus, has not led to new insights into the genetic determinants of 
Fhbl (Chouletefa/., 2010). 

Previous work has indicated that Fhbl plays a role in resistance 
against DON, a major virulence factor for F. graminearum. 
Lemmens etal. (2005) have shown that lines carrying the Fhb1 
resistant allele exhibit increased ability to conjugate DON into 
DON-3-glucoside (D3G). A highly DON-inducible barley uridine 
diphosphate-glycosyltransferase (UGT), capable of inactivating 
DON by transforming it to D3G, has been identified (Gardiner 
etal., 2010; Schweiger etal., 2010), and transgenic Arabidopsis 
carrying the barley UGT exhibits high DON resistance (Shin etal., 
2012). 

In contrast, a recent report has suggested that high levels of 
phenylpropanoids observed in F. gram/neari/m-inoculated spike- 
lets are responsible for the resistance conferred by Fhbl, which 
may be related to the also observed higher abundance of the 
biotic stress hormone jasmonic acid (JA; Gunnaiah ef a/., 201 2). In 
addition, Xiao etal. (2013) have identified JA signalling, but not 
biosynthesis, as specific for Fhb 7 -related resistance. Yet, this is 
difficult to reconcile with the results of studies using pure toxin 
(Horevaj etal., 2012; Lemmens etal., 2005). Possibly, in heavily 
inoculated spikelets, resistance is quickly overcome before the 
induction and function of UGT. Thus, induced UGT activity may be 
required for the inactivation of DON diffusion ahead of hyphae 
into the surrounding tissue. Alternatively, Zhuang etal. (2013) 
have associated the downregulation of a pectin methyl esterase 
inhibitor gene in strains lacking Fhbl with the QTL. Overexpres- 
sion in durum wheat of such a gene from Actinidia chinensis led to 
a reduction in FHB symptoms (Volpi etal., 2011). In addition, 
pathogenesis-related proteins and an ABC transporter have been 
found to be more abundant for the QTL (Gunnaiah etal., 2012), 
and an Rpg1 kinase homologue was mapped to the locus (Shen 
etal., 2006). 

This report investigates transcriptional differences of near- 
isogenic lines (NILs) segregating for Fhbl, and is the first to 
investigate Qfhs.ifaSA. We report a candidate gene - a lipid 
transfer protein (LTP) constitutively more abundant for the resist- 
ant Qfhs.ifaSA allele. A UGT gene more responsive in the Fhbl 
background - if Qfhs.ifaSA is also present - was identified and 
functionally tested by heterologous expression in yeast. 

RESULTS 

Resistance QTLs in a near-isogenic background confer 
different resistance levels 

In an effort to identify genetic determinants underlying two vali- 
dated FHB resistance QTLs in wheat, Fhbl and Qfhs.ifaSA, we 
generated NILs from a cross of the susceptible spring wheat cv. 



Remus and the FHB resistant line CM-82036, carrying the Fhbl 
and Qfhs.ifaSA QTLs. Four NILs from a BC 5 F 2 population were 
selected for the presence of the resistant alleles at both QTLs 
(NIL1), harbouring only the Fhbl (NIL2) or only the Qfhs.ifaSA 
(NIL3) resistance allele, or carrying the susceptible alleles at both 
QTLs (NIL4) (Fig. 1A). Diversity Arrays Technology (DArT) marker 
analysis showed that the percentage of Remus alleles varied 
between 96.3% (NIL1) and 99.8% (NIL4). 

To phenotypically assess the resistance levels, 150 heads per 
genotype were used to estimate disease progression, 12, 18 and 
24 days after point inoculation of single spikelets with F. gramine- 
arum spores. NIL4 developed a similar phenotype to the suscep- 
tible cv. Remus, and NIL1 exhibited the highest level of resistance 
among NILs, but did not perform as well as the parent CM-82036, 
which developed, on average, only one diseased spikelet (not 
shown). Intermediate resistance levels were observed in NIL2 
(Fhbi) and NIL3 (Qfhs.ifaSA) (f-test, P < 0.05; Fig. 1 B). At matu- 
rity, kernels of whole wheat ears were investigated for DON and 
D3G content, which confirmed the phenotyping results. The more 
resistant line NIL1 exhibited the lowest DON content (on average 
9.45 mg/kg) and the highest D3G/DON ratio (17%); NIL2 and NIL3 
contained, on average, 20.25 mg/kg and 18.52 mg/kg DON, 
respectively, and ratios of 12% and 14% D3G/DON. The suscepti- 
ble NIL4 showed the highest DON content (41 .35 mg/kg) and the 
lowest D3G/DON ratio (8%). 

However, inoculated samples from the time points used for the 
transcriptome analysis [72 h after inoculation (hai) and earlier] 
contained only trace amounts of DON. 

Transcriptional characterization of the wheat response 
to F. graminearum 

We assessed QTL-specific differential transcript abundance 
between the four NILs related to treatment (F. graminearum vs. 
mock) and time point (8, 24 and 72 hai) in three individual repli- 
cates using the Affymetrix 61 K wheat GeneChip. Several exam- 
ples are shown in Fig. 1C. In total, we identified 806 transcripts 
responding to F. graminearum infection in at least one time point 
when calculating the average (i.e. genotype-independent) differ- 
ences in transcript abundance between F. graminearum- and 
mock-inoculated samples (adjusted P < 0.05). The vast majority 
of these transcripts (792) were detected at 72 hai, less at 24 hai 
(62) and a few at 8 hai (29). The Fhbi- and £?ff?s./'fe-54-specific 
responses were generated by subtracting the F. graminearum! 
mock contrast of NIL4 (susceptible QTL alleles) from either NIL2 
(resistant Fhbi allele) or NIL3 (resistant Qfhs.ifaSA allele) 
(Fig. 1 A). The majority of transcripts identified with these contrasts 
showed a negative difference in response to F. graminearum, 
indicating a higher transcript accumulation in the presence of 
the susceptible QTL alleles (Table S1 , see Supporting Information). 
We confirmed the observed differences by reverse transcription- 
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Fig. 1 Characterization of near-isogenic lines (NILs) generated in this study. 

(A) Each NIL harbours a different composition of resistant (AA, BB) or suscep- 
tible (aa, bb) alleles of Fhb 1 and Qfhs.ifaSA. The Fusarium graminearum- 
specific responses for Fhbl and Qfhs.ifaSA were calculated by comparing the 
F. graminearumlmock contrasts of the indicated lines. Differences in constitu- 
tive gene expression were calculated by comparing mock-inoculated samples 
of the same lines. The 'interaction' contrast indicates nonadditive effects 
between Fhbl and Qfhs.ifaSA in NIL1 harbouring both resistance alleles. 

(B) Comparison of disease progression observed 18 days after inoculation in 
each of the NILs. (C) Heat map showing differential accumulation of 
exemplary transcripts between NILs at 72 h after inoculation (hai). The top 
panel represents Fusan'um-responsive transcripts: Ta. 12887. 1.S1_at, 

Ta. 23272. 1.S1_at and Ta.22732. 1 .S 1_at (all significantly changed for 
Qfhs.ifaSA). The bottom panel represents transcripts with constitutive 
differences in abundance: TaAffx.92377.1 .S1_at, Ta.6440.1.A1_at, 
Ta.1282.4.S1„at (all Qfhs.ifaSA) and Ta. 10270. 1 .A1_at (Fhbl). 

quantitative polymerase chain reaction (RT-qPCR) for a repre- 
sentative group of probe sets showing the highest relative 
changes. The investigated transcripts comprise F. graminearum- 
induced and QTL-specific genes (Fig. S1 A, see Supporting Informa- 



tion), and genes which are constitutively changed between mock- 
treated lines for one of the QTLs (Fig. S1B). Potentially, observed 
transcript levels can be confounded by single-feature polymor- 
phisms (SFPs) that lead to different binding affinities to the probe 
sets. However, such SFPs were detected for only a few transcripts. 
Transcripts were functionally annotated and compared against 
available databases, and with annotations reported previously by 
Jia ef al. (2009). For ease of comparison with Jia ef al. (2009), we 
classified the annotated genes by their putative roles as related to 
defence. In total, we identified 207 genes with a proposed role in 
defence and 264 involved in other categories (Table S2, see Sup- 
porting Information). 

Broad and late defence mechanisms are mounted in 
the absence of the Qfhs.ifaSA QTL 

The contrast to detect Qfhs.ifaSA gene-specific expression was 
calculated by subtracting the F. graminearumlmock contrast of 
NIL4 from the F. graminearumlmock contrast of NIL3 with the 
resistant Qfhs.ifaSA allele, and detected 352 transcripts as differ- 
ently regulated. Of these, 339 were detected at 72 hai and, with 
one exception (Ta.9975.2.S1_at, unknown), all showed a negative 
difference in response, indicating a higher induction or more sus- 
tained expression in the more susceptible genotype lacking Qfhs- 
ifaSA resistance. These 338 genes (Table S3, see Supporting 
Information) encode a variety of defence-related proteins. 

Twenty-one probe sets correspond to genes involved in the 
phenylpropanoid pathway, including multiple transcripts encoding 
phenylalanine ammonia lyases (PALs), the first dedicated step in 
the pathway. Only two more genes in this pathway that were not 
Qfhs.ifaSA specific were identified when calculating the average 
differences in transcript abundance between F. graminearum- and 
mock-inoculated samples across all lines. Most of the identified 
transcripts involved in the production of ethylene (ET) and JA 
showed higher accumulation for the susceptible Qfhs.ifaSA allele 
(one and nine transcripts, respectively). Similarly, F. graminearum- 
responsive transcripts potentially involved in DON detoxification 
were found to be predominantly more abundant in lines harbour- 
ing the susceptible Qfhs.ifaSA allele. These transcripts correspond 
to 12 UGTs (from a total of 14 UGTs detected), three MATE-type 
transporter genes (three) and nine pleiotropic drug resistance 
transporter genes (18). 

In contrast, transcripts corresponding to pathogenesis-related 
proteins were, in most cases, detected as a general response to the 
pathogen and not by the C?r7?s./ra-5.4-specific contrast. These 
included the subclasses of chitinases, glucanases, thaumatin and 
thaumatin-like proteins, as well as LTPs, which are involved in the 
early interaction with the pathogen (Walter ef al., 2010). However, 
the expression of germin-like proteins and protease inhibitors was 
almost exclusively detected in the more vulnerable line carrying 
the susceptible Qfhs.ifaSA allele (see also Table S2). 
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Table 1 Transcripts exhibiting a differential response for Qfhs.ifaSA* (adjusted P < 0.05, >2 fold change, FC) after inoculation with Fusarium graminearum 
spores. 

Significant 









FC, 
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FC, 


in other 


Mapped 


Affymetrix ID 


Annotationt 


Classification 


8hait 


24 hai 


72 hai 


contrasts§ 


toll 


Ta.9975.2.S1_at 


Unknown 


Unknown 


1.51 


0.95 


6.28 






Ta.14573.1.S1_at 


Hsp20/a-crystallin family protein (Oryza sativa) 


Defence 


0.23 


0.80 


1.68 






Ta.1 3859.1. A1_at 


OsCML16; calmodulin-related calcium sensor 
protein (Oryza sativa) 


Regulatory 


2.94 


0.64 


1.30 


-Fhb1 (24 hai) 




Ta.1 6827.1. S1_s_at 


Unknown 


Unknown 


0.75 


0.35 


1.25 






Ta.10311.1.S1_at 


Stress-responsive A/B barrel domain-containing 

nrntpin (Drvrz} 

UIULCMI \Uiyt-G juLIVu/ 


Unknown 


2.65 


1.01 


1.16 






Ta.1 3426.1. S1_s_at 


Os04g01 19500 protein (Oryza sativa) 
(UniRef90_Q7XTJ2) 


Unknown 


3.14 


0.89 


1.09 






TaAffx.64985.1.S1_at 


Putative uncharacterized protein (lea mays) 
(UniRef90_B4FQP2) 


Regulatory 


0.17 


2.02 


1.00 






Ta.9389.1.S1_at 


Seed maturation protein PM41 (Oryza sativa) 


Development 


2.20 


0.31 


0.98 






Ta.20205.1.A1_s_at 


VER2 (Triticum aestivum) (UniRef90_O80370) 


Development 


1.14 


0.48 


0.95 






Ta.24801.1.S1_at 


Hsp20/a-crystallin family protein (Oryza sativa) 


Defence 


0.25 


0.95 


0.90 






Ta.5004.1.S1_at 


Thiol protease SEN 102 precursor (Oryza sativa) 


Metabolism 


0.26 


1.70 


0.89 


-Fhbl (8 hai) 




Ta.5004.3.S1_a_at 


Thiol protease SEN 102 precursor (Oryza sativa) 


Metabolism 


0.47 


1.19 


0.88 


-Fhbl (8 hai) 




TaAffx.1 12045.1. 51_x_at 


Glutathione-5-transferase (Oryza sativa) 


Defence 


3.20 


1.23 


0.03 


+Fhb1 (8 hai) 




Ta.5610.1.S1_at 


Nonsymbiotic haemoglobin 2 (Oryza sativa) 


Miscellaneous 


7.01 


0.40 


0.02 




4BL, 4DL 



*The expression ratio was calculated from the F. graminearumlmock response observed in NIL3 (resistant Qfhs.ifaSA allele) and the response provided by NIL4 
(susceptible allele). Genes more abundant in the absence of Qfhs.ifaSA at 72 hai are excluded here and provided in Table S3 (see Supporting Information). 
tAnnotations were assigned if the £ value from the BlastX results was r£l0~ 10 ; otherwise the transcripts were denominated as 'unknown'. 
^Differences in response are the means of three biological replicates. Statistically significant values are given in bold type. 
§-/+ indicates significant positive/negative differences in response for the indicated contrast. 

IMap locations were identified by searching the mapped expressed seguence tag (EST) database of Graingenes for the probe set seguences using BlastN and an E cut-off 
of <10- 10 . 

At time points other than 72 hai, five probe sets with a positive 
difference in response to F. graminearum were identified for 
Qfhs.ifaSA (8 hai), indicating higher accumulation for the resist- 
ant allele. These encoded a stress-responsive A/B barrel domain- 
containing protein, a calmodulin-related calcium sensor protein, a 
glutathione-5-transferase (GST), a nonsymbiotic haemoglobin and 
one unknown protein. Notably, two of these genes also exhibited 
a negative difference at 72 hai, with higher accumulation in the 
genotype with the susceptible Qfhs.ifaSA allele (Ta. 561 0.1 .S1_at, 
nonsymbiotic haemoglobin 2 and TaAffx.1 1 2045.1 .S1 _x_at, GST). 
Fusarium graminearum-responswe transcripts exhibiting a signifi- 
cant difference in response for Qfhs.ifaSA are summarized in 
Table 1 . 

An LTP is constitutively more abundant in the 
presence of Qfhs.ifaSA 

Twenty-eight genes were constitutively more and 63 less abundant 
in the presence of the resistant allele when comparing mock- 
inoculated samples of NIL3 and NIL4 differing in Qfhs.ifaSA (Table 
S4, see Supporting Information, and Table 2, which summarizes the 
group of more abundant transcripts). Of these, 16 probe sets 
corresponded to defence-related genes, including three heat shock 
proteins more abundant only at 8 hai (two isogenes) or 72 hai (one 
isogene). The largest difference in transcript abundance was 
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observed for Ta.1 282.4. S1_at, which corresponds to an LTP. The 
probe set derives from an expressed sequence tag (EST) (GenBank 
accession CA635994), which is identical to the gene model of 
TaLTP9.4a (GenBank accession AJ852541, cv. Apache; Boutrot 
era/., 2007). We used specific primers for the 5' untranslated 
region (UTR) published forTaLTP9.4a to map Ta.1 282.4.S1 _at to 
the long arm of chromosome 5A using a collection of Chinese 
Spring nulli-tetra lines. Deletion bin mapping of this gene using 
four Chinese Spring deletion lines for 5AS and five deletion lines for 
5AL produced the specific PCR fragment for all lines investigated, 
suggesting that this gene is located on the pericentromeric region 
of 5AL. The gene model of this LTP contains one 440-bp intron, 
which separates the 369-bp coding sequence (CDS) into two exons 
of 360 and 9 bp. Resequencing of cloned PCR products from both 
the resistant CM-82036 and the susceptible Remus identified two 
highly similar LTP genes, which exist in both genotypes. These 
genes, potentially paralogues positioned on the same locus, differ 
in two and six single-nucleotide polymorphisms (SNPs), respec- 
tively, from the EST template (Fig. S2, see Supporting Information). 

Fusarium graminearum-responsive 
Fhb 1 -specific transcripts 

The F. graminearum-spectik and Fhb 1 -dependent response was 
generated by subtracting the F. graminearum/mock contrast of 
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Table 2 Transcripts constitutively more abundant (adjusted P £ 0.05, >2 fold change, FC) in the presence of the resistance alleles of either Fhbl or QfhsJfaSA*. 
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Unknown 


Unknown 


30.91 


21.96 


36.91 


5 


Ta. 13161.1 .S1_at 


Putative extensin-like cell wall protein (Jia et aL, 2009) 


Cell wall 


21.24 


23.13 


36.24 


8,9 


Ta.10270.2.A1_a_at 


Wax synthase (Oryza sativa) 


Metabolism 


5.62 


6.42 


6.1 
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Putative uncharacterized protein Sb02g001860 


Unknown 


5.15 


5.18 


4.94 
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2.89 


2.65 


2.99 
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Unknown 


Unknown 


2.54 


2.41 


2.84 




Ta.5633.1.S1_a_at 


Unknown 


Unknown 


1.02 


3.76 


2.7 




Ta.25791.1.A1_at 


Unknown 


Unknown 


3.15 


4.83 


2.55 


10 
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Putative uncharacterized protein (Oryza sativa Indica 


Unknown 


1.86 
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1.89 
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Unknown 
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1.39 


2.25 


1.39 
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Seed maturation protein PM41 (Oryza sativa) 


Development 


0.33 


2.67 
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TaAffx.616.2.S1_s_at 


Hypothetical protein (Jia etal., 2009) 


Unknown 


0.58 


2.83 
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Subtilisin— chymotrypsin inhibitor 2 (Jia etal., 2009) 
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Expressed protein (Oryza sativa) 


Unknown 


9.92 


19.24 


17.4 
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RCN1 Centroradialis-likel homologous to TFL1 gene; 


Regulatory 


3.2 
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contains Pfam profile PF01 161: 
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Ta.1 1043.1.A1_at 


Unknown 


Unknown 


5.28 


6.67 


5.8 




TaAffx.70660.1.S1_at 


Hsp20/a-crystallin family protein (Oryza sativa) 


Defence 


1.5 


1.03 


4.77 




TaAffx.92377.1.S1_at 


Predicted protein (Populus trichocarpa) 


Regulatory 


3.74 


4.23 


4.41 






(UniRef90_B9GJ95) 












TaAffx.9181.1.S1_at 


Unknown 


Unknown 


3.61 


3.13 


4.19 


7 


Ta.5776.1.S1_at 


Unknown 


Unknown 


0.96 


1.03 


3.81 




Ta.25955.1.A1_at 


Unknown 


Unknown 


2.7 


1.93 


2.82 




Ta.23102.1.S1_at 


HEAT repeat family protein (Oryza sativa) 


Miscellaneous 


5.4 


2.37 


2.64 


2DS, 5 AS, 


TaAffx.23477.1.S1_at 


Unknown 


Unknown 


2.53 


2.62 


2.53 




Ta.6352.3.A1_a_at 


XH domain-containing protein (Oryza sativa) 


Regulatory 


2.42 


2.7 


2.51 




TaAffx.86453.1.S1_s_at 


Unknown 


Unknown 


1.92 


2.89 


2.24 




Ta.5597.1.S1_at 


Auxin-repressed protein (Oryza sativa) 


Unknown 


1.62 


1.72 


2.15 


7 4DL, 6BS 


Ta.12038.1.A1_at 


MDR-like ABC transporter (Oryza sativa) 


Putative trichothecene 


1.86 


2.13 


2.06 








detoxification and transport 










Ta.5281.1.S1_at 


Hydrolase, a/fj fold family domain-containing protein 


Metabolism 


3.42 


3.91 


1.85 






(Oryza sativa) 












Ta.4381.2.S1_a_at 


Unknown 


Unknown 


2.5 


2.13 


1.74 


7AS 


Ta.29578.2.S1_at 


Nonsymbiotic haemoglobin 2 (Oryza sativa) 


Miscellaneous 


0.84 


2.06 


1.7 




Ta.867.1.A1_at 


Nucleotide pyrophosphatase/phosphodiesterase 


Metabolism 


2.39 


2.41 


1.68 






(Andropogoneae) (UniRef90_B6U857) 












Ta.13156.1.S1_at 


Unknown 


Unknown 


2.13 


1.96 


1.64 




Ta.8274.1.A1_at 


Unknown 


Unknown 


2.37 


2.26 


1.56 




Ta.540.1.Sl_at 


Unknown 


Unknown 


2.02 


1.34 


1.53 


1AL, 3 AS 


Ta.19896.1.A1_at 


Os09g0425200 protein (Oryza sativa) 


Unknown 


2.24 


1.48 


1.36 






(UniRef90_Q69PC6) 












Ta.24801.1.S1_at 


Hsp20/a-crystallin family protein (Oryza sativa) 


Defence 


4.55 


1.03 


1.36 




Ta.1 1 1 1 1 .1 .A1_at 


RCN4 Centroradialis-likel homologous to TFL1 gene 


Development 


2.92 


1.81 


1 






(Oryza sativa) 












Ta.14573.1.S1_at 


Hsp20/a-crystallin family protein (Oryza sativa) 


Defence 


4.77 


1.29 


0.86 




*The expression ratios were calculated by comparing mock-treated samples of NIL2 (resistant Fhbl allele) or NIL3 (resistant QfhsJfaSA allele) with NIL4 (susceptible alleles). 
tAnnotations were assigned if E from the BlastX results was s10" w ; otherwise, the transcripts were denominated as 'unknown'. 
♦Differences in transcript abundance are the means of three biological replicates. Statistically significant values are given in bold type. 
§Numbers in the single-feature polymorphism (SFP) column indicate for which probes SFPs were detected. 

UMap locations were identified by searching the mapped expressed sequence tag (EST) database of Graingenes for the probe set sequences using BlastN and a E < lO" 10 . 
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NIL4 (both susceptible QTL alleles) from the F. graminearum/mock 
contrast derived from NIL2 (resistant Fhb1 allele). We identified 
16 transcripts exhibiting a differential response for Fhb 7 (Table 3). 
At 8 hai, only one transcript exhibited a positive difference 
in response (TaAffx.1 1 2045.1 .S1_x_at; GST), and two probe 
sets had a negative difference in response (Ta.5004.1.S1_at; 
Ta. 5004.1 .S1 _x_at; thiol protease precursors). With one exception 
(Ta.14071.1.S1_x_at; integral membrane protein), all of the tran- 
scripts identified at 72 hai (12) and 24 hai (two) showed a nega- 
tive difference in response, indicating less accumulation in the line 
carrying the resistant Fhb1 allele compared with the line carrying 
the susceptible allele. 

Fifty-four transcripts exhibited constitutive Fhb 1 -dependent dif- 
ferences in abundance, when comparing mock-inoculated samples 
of NIL2 and NIL4. Most transcripts showed a consistent expression 
pattern throughout all time points (Table S5, see Supporting Infor- 
mation). Of these 54 transcripts, 17 were constitutively more 
abundant in the presence of the resistant Fhb1 allele (Table 2). 
These transcripts encode, for instance, an extensin-like cell wall 
protein, a wax synthase, a calmodulin-related calcium sensor 
protein, a subtilisin-chymotrypsin inhibitor (only detected for 
24 hai), two seed maturation proteins and 1 1 proteins of unknown 
function. 

To identify Affymetrix probe sets mapping in close vicinity to the 
physical positions of markers linked to Fhb1, we employed the 
available partial genomic sequence of chromosome 3B of cv. 
Chinese Spring (Choulet etal., 2010), which included the region 
harbouring Fhbl.We located markers sts32, sts66, sts194, sts142, 
sts80, sts189 and umnW (Liu etal., 2006, 2008) on ctg0954b of 
3BS (GenBank accession FN564434, Fig. 2). The markers span a 
distance of 2.94 megabases (sts66 to sts189) and the region 
comprising the putative Fhb1 region spans about 274 kb (sts32 
to umnW). From the genes represented on the 61 K wheat 
GeneChip, we mapped 35 with at least 98% identity to ctg0954b. 
Four of the 35 mapped transcripts were changed significantly 
for Fhb1 in our experiments, including TaAffx.1 2498. 2. S1_at 
and Ta.22694.1.A1_at (constitutively more abundant for the 
resistant Fhb1 allele, both unknown) and Ta.28185.1.S1_at and 
Ta.6066.2.S1_a_at (constitutively more abundant for the suscep- 
tible allele, splicing factor and unknown, respectively; Fig. 2). The 
majority of transcripts changed for Fhb1 could not be mapped to 
the interval. Only eight GeneChip probe sets correspond to the 54 
annotated genes on the contig. 

Interaction effects between both QTLs identify 
additional Fhb 1 -associated candidates 

We also investigated whether the combined QTLs in NIL1 
exhibited interaction effects between Fhbl and Qfhs.ifaSA 
when compared with the transcript accumulation observed for 
single QTL contrasts. These were calculated by subtracting the 



F. graminearumlmock contrasts observed in NIL2 [Fhbl) and 
NIL3 (Qfhs.ifa-5A) from the contrast observed in NIL1 (both resist- 
ant alleles; interactions = NIL1 - NIL2 - NIL3 + NIL4; adjusted 
P < 0.05). The pathogen-responsive and constitutively changed 
transcripts are given in Tables S6 and S7, respectively (see Sup- 
porting Information). 

Twenty-three of the F. graminearum-respons\\/e transcripts that 
showed interaction effects at 72 hai were similarly abundant at 
high levels in both lines lacking the 5A resistance (NIL2 and NIL4). 
Hence, the interaction contrast is generated by the difference 
between NIL1 and NIL3, indicating potential Fhbl dependence (for 
instance, Ta.1 2887.1 .S1_at; Fig. 1 C). We tested whether these 
genes were also changed significantly for Fhbl in the genetic 
background with the resistant Qfhs.ifaSA allele by calculating the 
difference in response to F. graminearum between NIL1 and NIL3 
(which also includes potential interaction effects; Table S8, see 
Supporting Information). This contrast identified the same three 
UGTs, two GSTs, two P450 monooxygenases and a laccase 
protein, as well as three metabolism-related proteins and four 
unknown. The annotated gene with the highest difference in 
response was TaAffx.1 1 1 585.1 .S1 . A similar contrast was also 
generated for Qfhs.ifaSA in the background of Fhbl (NIL1 - NIL2; 
Table S9, see Supporting Information), which, however, did not 
yield additional candidates. 

The Fhb /-associated TaUGT12887 confers weak 
resistance to DON 

BlastN analysis identified TaAffx.1 11 585.1 .S1_at as well as 
Ta.1 2887.1. S1_at as close homologues to HvUGT13248, a barley 
gene encoding a UGT with a proven ability to inactivate DON 
(Schweiger et al., 2010). As resistance conferred by Fhbl has been 
shown previously to be correlated with greater ability to gluco- 
sylate DON (Lemmens etal., 2005), both probe sets were attrac- 
tive starting points for further investigations. Reconstruction of 
the gene models using the publically available genomic 454 
sequences (http://www.cerealsdb.uk.net) of cv. Chinese Spring 
showed that the ESTs corresponding to both probe sets originated 
from the same gene. The assembled gene model contained two 
introns, which separated a large first exon (705 bp) and two 
smaller exons of 385 and 338 bp in length. The encoded open 
reading frame (ORF) (1428 bp) was designated TaUGT12887 
(GenBank accession JX624788).We used a collection of nulli-tetra 
and ditelosomic lines to map TaUGT12887 to chromosome 5BL. 
Resequencing of PCR products from both parental genotypes of 
the mapping population showed that the coding sequence was 
identical in Chinese Spring and the two parental lines used in this 
study, CM-82036 and Remus. 

To test whether TaUGT12887 encodes a protein with the ability 
to inactivate DON, we recoded the sequence for optimal expres- 
sion in yeast. The toxin-sensitive Saccharomyces cerevisiae strain 



© 2013 THE AUTHORS. MOLECULAR PLANT PATHOLOGY PUBLISHED BY JOHN WILEY & SONS LTD AND BSPP 
MOLECULAR PLANT PATHOLOGY (2013) 14(8), 772-785 



778 W. SCHWEIGER ef al. 



03 £J 
k_l s/1 



rsi 



O 
rsi 

o 



CO ' — ' — 



^ E 

O „ =3 

V 1? 
so -e 03 

3 "> -a 
° to S 

itv 

o--a ~ 

£ S-s 

ru OJ O 
_Q Q_ _E 
p X H3 

is; 

OJ o_ k_j 
iE O 





















_j 




_i 


_J- 

Q 


















^ . 






rn 


LO 












































CO 


CO 




co 


CO 


















*— 






m 


on 




















r 
























< 


< 




< 


CQ 
























m 






'rO 


ra 


re 


ra 




ra 


CD 




'ro 


'ro 


'ro 


'ra 


oo 






CO 


oo 


OO 


oo 




oo 




oo 


rsi 


oo 


? 




























? 


? 


? 


? 


? 


? 


? 


? 


? 


? 


? 




42 




42 


42 


42 


42 


42 


42 


42 


42 


42 


42 


42 


i*i 


























tg 


eg 


eg 


eg 


jg 


i 


eg 


eg 


eg 


eg 


eg 


jg 


jg 


o 


O- O- Cr O- O 




Q 


O 


Q 





rsi 1-0 ^ 

rn ^ ^ ^ 



ci ^ 



s £ 6 S 



i ^ 



*-> r-s 



cn 



^- <— H o ^" 

sT st ^ lo in i- 



sO rsj 
sO 

>< >< 

< < 



o o o o o o 



^ ^ ^ o o o ^- ^- ^- ^ ^ ^~ 



© © o ^ o o 



SEE a E 

_Q_Q_Q y y.>+=i_Q 

SJSiSojsu'S-giS 

cu cu <u ^ ^ ±± ^ OJ 

5 5 5 q q o_ 5 



o -g 



O O o <2 



2 S n: 

d □ *s .E 

t_j sj sb aj 

£H £rj fu o 

rsi rs! ^ 

5m 9 i 



2 p 

'I & 

03 

i~ ro 



5 5 5 



CO .-ts £= 



dj O 



o o o a> .2 ■— 

(U U ,^ 

o o "a t± m 



tu ra P 

Kj iZ <o OJ - 

- s. 



O aj 03 

!S |.-s3 



<£s5tl !£ 



'st 



i i : 



ix> r- t— 

if co in co 

^ CO rsi <- 

(T3 iT3 fD 03 



m ' — 
CO <— 



Ln m rsi 
r- m 
Ln ^ m 



S 5 



o 



cn 



=5 



^ vi ^ s 

C Q_ O 

S iJJ Ci Q_ 

o cu s— on 

5 3 -= 



m- -£= "O 
— ^ ~ OJ 



3 -Si 



ST IS 5 -a 



MOLECULAR PLANT PATHOLOGY (2013) 14(8), 772-785 
2013 THE AUTHORS. MOLECULAR PLANT PATHOLOGY PUBLISHED BY JOHN WILEY & SONS LTD AND BSPP 



RNA profiling of Fusarium resistance QTLs 779 



probe set ID 




SSR marker 



sts66 

0.9 cM, 268 kb 



TaUGT12887 HvUGT13248 



vector control 



sts194 



0.4 cM, 1055 kb 



sts142 



1.2 cM, 1571 kb 



TaAffx.12498.2.S1_at 
Ta.22694.1.A1_at 

Ta.6066.2.S1_a_at 



Ta.28185.1.S1_at 



sts32 



OcM, 2193 kb 



sts80, umn10 

0.4 cM, 2461, 2467 kb 




sts189 



2945 kb 



Fig. 2 Comparison of transcripts with the physical map of the Fhb 1 locus. 
Simple sequence repeat (SSR) markers mapped to the sequence of the 
Chinese Spring BAC library contig ctg0954b of chromosome 3BS are given 
with genetic distances to each other and their physical position. White bars 
represent positions of annotated genes. Black bars represent mapped 
expressed sequence tag (EST) sequences corresponding to the Affymetrix 
wheat GeneChip. Probe set IDs of transcripts differentially changed for Fhb1 
are indicated. 

YZGA51 5 was transformed with a plasmid containing the custom 
synthesized ORF under the control of the strong constitutive ADFil 
promoter. Spotting assays showed that TaUGT12887 conferred 
only weak resistance against DON compared with the barley 
DON-inactivating HvllGT13248 (Fig. 3), although the proteins 
were highly expressed in cell extracts, as confirmed by immuno- 
detection of the N-terminal c-Myc tag (Fig. S3, see Supporting 
Information). 
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Fig. 3 Functional evaluation of the fM> /-responsive uridine 
diphosphate-glycosyltransferase (UGT) TaUGT12887. Spottings of 
toxin-sensitive yeast strains expressing TaUGT12887, HvUGT13248 as a 
positive control or the empty plasmid as negative control. Two individual 
transformants of each strain were spotted in different dilutions on medium 
containing increasing amounts of toxin. DON, deoxynivalenol. 

DISCUSSION 

We investigated the effect of the Qfhs.ifaSA QTL on gene expres- 
sion and reinvestigated the host response conferred by Fhb1 using 
NILs harbouring either, both or none of the resistance alleles of 
Fhb1 and Qfhs.ifaSA. The NILs have the European spring wheat 
cultivar Remus as the recurrent parent. Previous reports have 
shown that the impact of both QTLs is highly dependent on the 
genetic background of the recurrent parent (von der Ohe etal., 
2010; Salameh etal., 2011). In the present study, we observed 
that Qfhs.ifaSA (type I resistance) reduces disease severity at a 
similar level to lines with Fhb1 (type II resistance) after point 
inoculation — a method more sensitive to the assessment of type II 
resistance. The high level of resistance conferred by Qfhs.ifaSA 
may have been more evident with the cautious single floret inocu- 
lation method used. Most of the F. graminearum-responswe genes 
found represent a general response to the pathogen, regardless of 
resistance level and QTL. 

Previous transcriptomic studies on the F. graminearum-wheat 
interaction have led to the characterization of multiple defence- 
associated gene classes that describe a general or genotype- 
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specific response of wheat to F. graminearum or DON (Ansari 
etal., 2007; Bernardo etal., 2007; Diethelm etal., 2012; 
Ding etal., 2011; Foroud etal., 2011; Golkari etal., 2009; 
Gottwald etal., 2012; Steiner etal., 2009; Walter etal., 
2008). Walter ef al. (2008) examined the gene expression of DON- 
treated doubled haploid (DH) lines segregating for Fhbl, and 
reported the Fhbl -dependent upregulation of several genes. Jia 
etal. (2009) investigated a NIL pair differential for the resistant/ 
susceptible allele of Fhbl. They found 14 transcripts exhibit- 
ing a differential expression after comparing F. graminearum- 
inoculated lines directly. Of these 1 4 transcripts, we identified five 
in our study: Ta.23324.1 .S1_at (unknown), Ta. 22694.1. A1_at 
(unknown), Ta.28185.1 .S1_at (splicing factor), Ta.438.3.S1_at 
(nuclease) and Ta.9489.1 .S1 _at (regulatory gene). In addition, only 
a few genes identified here as QTL specific were identified in other 
comparable studies also using the Affymetrix GeneChip. The 
F. graminearum-responsH/e Ta. 2471 5.1. S1_at, encoding a peroxi- 
dase, was more abundant for Fhbl at the earliest observed time 
point of 8 hai and was also found to be upregulated in response to 
DON, but only for the susceptible genotype Superb (24 hai; Foroud 
ef al., 2011). This peroxidase might play a role in the generation of 
reactive oxygen species (ROS). Ding etal. (2011) also observed 
earlier abundance of peroxidases in the resistant Wangshuibai 
compared with a susceptible mutant, and speculated that the 
production of ROS induces JA and ET biosynthesis (which was not 
observed for Fhbl in this work). In addition, Xiao etal. (2013) 
reported JA signalling to be upregulated for the QTL. 

The complex data resulting from transcriptomic studies may 
become clearer if potential biochemical mechanisms are also con- 
sidered. A reasonable basis for resistance to initial infection could 
be a result of a preformed barrier or antifungal protein or plant 
metabolite. From the candidates listed in Table 2, the 50-70-fold 
more highly expressed LTP is the obvious candidate gene. With 
regard to Fhbl, a strong case has been made that the spread of 
resistance is caused by resistance to the toxin DON, which is 
necessary for the spread of Fusarium. 

Potential links to FAM-regulated DON detoxification 

Although, in the highly resistant CM-82036, the spread of the 
fungus is impeded in the initially inoculated spikelet, in the NILs, 
the fungus manages to overwhelm resistance responses and infect 
the rachis, where cell wall appositions create a barrier against 
spread to the next spikelet. The synthesis of phenylpropanoids 
seems to play an important role (Gunnaiah etal., 2012), but the 
fungus seems able to inhibit its formation by the production of the 
protein biosynthesis inhibitor DON (Jansen etal., 2005). DON 
detoxification is not necessarily constitutive, as multiple UGTs 
have been shown to be inducible by DON (Gardiner etal., 2010; 
Poppenberger ef al., 2003; Schweiger ef al., 201 3). Thus, rapid and 
high inducibility of such UGTs, before the toxin can block protein 



synthesis, may be required for resistance. Even NIL1 had, on 
average, four infected spikelets, indicating that the fungus can 
overwhelm resistance mechanisms in the inoculated spikelet. Pre- 
sumably, in the more susceptible lines NIL2 and NIL4, the resist- 
ance mechanisms in the sampled inoculated spikelets have 
already been overwhelmed and QTL-specific differences are no 
longer evident. We found one candidate UGT, designated 
TaUGT12887, in the Fhbl background, but only if basal resistance, 
caused by the presence of Qfhs.ifaSA, is already high (Fig. 1 C). 
Potentially, interaction effects between both QTLs could also play 
a role in the greater abundance of this gene. Yet, when expressed 
in toxin-sensitive yeast, TaUGT12887 only exhibits weak activity 
against DON compared with HvUGTI 3248 from barley (Schweiger 
etal., 2010). 

In Brachypodium distachyon, all members of the orthologous 
gene cluster to TaUGTI2887 were found to be induced by DON, 
but only two of six of these genes coded for the ability to inacti- 
vate the toxin (Schweiger etal., 2013). A similar situation can be 
expected for the orthologous gene families on the three genomes 
of wheat, wherein TaUGT12887 might be a DON-responsive, but 
weakly active UGT, whereas yet uncharacterized paralogues, 
potentially not represented on the GeneChip, perform the primary 
detoxification function. 

The formation of glutathione adducts may also contribute to 
DON detoxification. One GST in our data exhibited a positive 
difference in response to F. graminearum inoculation for Fhbl 
(8 hai, TaAffx.11 2045.1 .S1_x_at): This probe set corresponds 
to a putative tau-class GST, which was also identified by 
other contrasts. Two further GSTs (TaAffx.1 2851 4.1 .S1_at and 
Ta.886.2.S1_at, 72 hai) were identified in the additional Fhbl 
contrast, in the background of Qfhs.ifaSA. Transcriptomic studies 
in DON-infected barley indicate that the F. graminearum virulence 
factor DON most probably depletes the glutathione pool. The 
involvement of GSTs in the detoxification of DON is still unproven, 
but the nonenzymatic formation of DON-glutathione adducts has 
been shown (Gardiner etal., 2010). Recently, evidence for the 
occurrence in planta of a DON-glutathione conjugate and derived 
processing products in wheat has been demonstrated (Kluger 
etal., 2012). 

Lemmens etal. (2005) have proposed that the Fhbl locus 
either encodes a UGT or a regulatory gene affecting detoxification 
capability, but neither the UGT nor the GST found to be differen- 
tially expressed in this study is located in the Fhbl interval. With 
regard to a potential regulator, a gene of unknown function, 
Ta. 22694. 1.A1_at, is interesting. It is constitutively more abun- 
dant (about 30-fold) at all observed time points in the line har- 
bouring the resistant Fhbl allele. This gene is located close to the 
predicted location of Fhbl between sts!42 and sts32 (see Fig. 2). 
No homologous DNA is present in Sorghum bicolor, rice or B. dis- 
tachyon. Only three ESTs corresponding to this sequence, two 
sense transcripts (cv. Chinese Spring and cv. Stephens) — for which 
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the probe sets have been designed — and one antisense transcript 
(cv. Chinese Spring), are present in the database. It is unclear 
whether there is a meaningful (short) ORF or whether the tran- 
script has a regulatory role. 

A constitutively more highly expressed LTP associated 
with Qfhs.ifaSA 

We identified 338 genes exhibiting a negative difference in 
response for F. graminearum at 72 hai for Qfhs.ifaSA, indicating 
higher accumulation in the presence of the susceptible Qfhs.ifaSA 
allele. One hypothesis that explains the first counterintuitive obser- 
vation of a higher upregulation of multiple defence-related genes 
in the susceptible cultivar is that the fungal biomass is higher in the 
lines lacking the Qfhs.ifaSA resistance allele. Consequently, more 
plant tissue is undergoing a defence response at a given time. 

Of the few candidate genes for Qfhs.ifaSA proposed in the 
literature, none has been identified in this study. Steiner etal. 
(2009), using a cDNA-amplified fragment length polymorphism 
(AFLP) approach, investigated the QTLs in lines harbouring both 
Fhb1 and Qfhs.ifaSA, and identified transcripts encoding a UGT 
and a PAL upregulated at 6 hai. Neither of these genes has been 
identified as specific for either QTL in our study. 

A promising Qfhs.ifaSA candidate gene is Ta.1282.4.S1_at, a 
type 1 nonspecific LTP. This gene is constitutively more abundant in 
lines with the Qfhs.ifaSA resistant allele, exhibiting one of the 
largest differences in abundance in the microarray data. The posi- 
tion of Ta. 1282.4. S1_at at 5AL, in close proximity to the centro- 
mere, is in the current confidence interval for Qfhs.ifaSA 
(Buerstmayr etal., 2003). Foroud etal. (2011) also reported 
Ta.1 282.4. S1_at to be constitutively more abundant in a DH line 
with CYMMIT1 1, a FHB resistance donor with an unknown resist- 
ance QTL in its pedigree. However, screening for Qfhs.ifaSA using 
the flanking markers Gwm293 and Gwm156 did not identify 
Qfhs.ifaSA in their studied genotypes. The same LTP gene identi- 
fied in this study was found to be associated with Qfhs.ifaSA in a 
microarray study comparing lines from a cross of Wuhanl and 
NuyBay (T. Ouellet, Agriculture and Agri-Food Canada, Ottawa, 
personal communication). 

Multiple roles in defence against pathogens are ascribed to LTPs 
(Yeats and Rose, 2008). Sun etal. (2008) demonstrated that a 
different recombinant LTP protein from wheat can inhibit the 
in vitro growth of multiple fungal wheat pathogens, including 
F. graminearum. The Arabidopsis thaliana DIRI gene encodes an 
LTP that binds to lipids and elicits long-distance signalling, and 
thus contributes to systemic acquired resistance (Maldonado 
etal., 2002). Another member of the large wheat LTP family has 
also been shown to compete for plasma membrane receptors with 
elicitins, proteins secreted by fungal pathogens from the Phytoph- 
thora or Pythium genera that elicit hypersensitive responses 
(Buhot etal., 2001). The early interaction of F. graminearum and 



the host is characterized by the release of lytic enzymes which 
penetrate directly the outer epidermal layers (Kang and 
Buchenauer, 2002). The constitutively higher abundance of LTP, 
also reported to act in cutin deposition (Yeats etal., 2010), might 
coincide with the type I resistance against initial fungal penetra- 
tion attributed to Qfhs.ifaSA. LTPs might also contribute to toxin 
resistance, as the overexpression of an A. thaliana LTP leads to 
higher resistance to the trichothecene toxin trichothecin 
(McLaughlin etal., 2011). Further work to functionally character- 
ize this LTP by stable overexpression in wheat is planned. If the 
role of the LTP as an F. graminearum resistance determinant can 
indeed be confirmed, new questions related to food safety must be 
addressed. Many PR-14 class LTPs do not only have a role in plant 
defence, but are also relevant food allergens (Palacin ef al., 2007). 

Taken together, this study presents candidate genes associated 
with the F. graminearum resistance QTLs Fhbi and Qfhs.ifaSA 
and functionally validates an Fhbi -related UGT, which improves 
resistance to DON in toxin-sensitive yeast to a limited extent. With 
the release of the wheat genomic sequence in the foreseeable 
future, these candidates and available transcriptional data can be 
mapped and compared with the draft sequence of the (suscepti- 
ble) cultivar Chinese Spring. Positional cloning efforts in resistant 
germplasm, combined with functional validation of candidates in 
wheat using TILLING (targeting induced local lesions in genomes) 
populations or virus-induced gene silencing, will provide further 
insights into the resistance mediated by both QTLs. For Qfhs.ifa- 
SA, additional efforts are required to resolve the locus in the 
centromeric region of chromosome 5A. 

EXPERIMENTAL PROCEDURES 
Development of NILs 

NILs were generated via five backcrosses to cv. Remus (FHB susceptible, 
Sappo/Mex//Famos) of a DH line identified from a CM-82036 (FHB resist- 
ant, Sumai3/Thornbird) by Remus population. Selfed F 2 progeny of BC 5 
plants were screened for the presence of the homozygous resistance QTL. 
NIL1 possessing the resistance alleles of Fhbi and Qfhs.ifaSA, NIL2 pos- 
sessing Fhbi, NIL3 carrying the Qfhs.ifaSA resistance allele and NIL4 
carrying the susceptible alleles at Fhbi and Qfhs.ifaSA were selected 
(Fig. 1A).The presence of the resistance QTL allele on 3BS was verified by 
the diagnostic co-segregating marker umnW.lhe genotype of Qfhs.ifaSA 
was verified with flanking simple sequence repeat (SSR) markers barc186, 
gwm1057, barc56, gwm293, bard 17 and gwml29 on the short arm of 
chromosome 5A, and bard, barcl80, barc40 and gmwl56 on the long 
arm of 5A. The four NILs, CM-82036 and Remus were also genotyped with 
DArT markers by Triticarte Pty. Ltd. (Canberra, Australia; Akbari etal., 
2006; Wenzl etal., 2004). 

Glasshouse experiments and tissue sampling 

The glasshouse experiment and inoculum production were conducted as 
described in Steiner ef al. (2009) with modifications. The experiment was 
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performed in three replications and the lines were challenged by inocula- 
tion with 10u.L of a 20 000 spores/mL F. graminearum suspension or 
water (about 200 conidia). Five heads per genotype, treatment and repli- 
cation were harvested at 8, 24 and 72 hai. Treated spikelets were sepa- 
rated into lemma, palea and the subtending section of the rachis and the 
reproductive tissues, immediately frozen in liquid nitrogen and stored at 
-80 °C. For phenotype scoring, one central spikelet (two florets) of 50 
wheat heads per genotype was inoculated on the field with 100 conidia 
per floret, and diseased spikelets were counted after 12, 18 and 24 days. 
The inoculum concentration used for phenotyping allowed better differ- 
entiation of phenotypes between NILs. Statistically significant differences 
between all lines were calculated using pairwise f-tests (adjusted 
P < 0.05). The experiment was repeated twice. 

RNA preparation and microarray hybridization 

Sampled lemma, palea and rachis from five heads per experimental con- 
dition were pooled and ground to a fine powder using a mixer mill (type 
MM 301, Retsch, Dusseldorf, Germany). RNA was extracted according to 
Steiner ef al. (2009). We employed the Affymetrix wheat Genome Array 
(Affymetrix, Santa Clara, CA, USA) for the identification of differential 
transcript abundances in our samples. The generation of labelled cRNA, 
hybridization and data acquisition were performed at the Microarray 
Facility of the University of Minnesota following the standard protocol 
described by Boddu era/. (2006). 

Array and SFP analysis 

Statistical analysis was performed using the R statistical analysis environ- 
ment (http://www.r-project.org) and packages of the Bioconductor suite 
(http://www.bioconductor.org). We used the 'gcrma' package (Wu era/., 
2004), 'vsn' (Huber ef al., 2002), the fitPLM function ('affyPLM'; Bolstad, 
2004) and the ImFIT function of 'limma' (Smyth, 2004). To identify signifi- 
cant differentially expressed probe sets, we calculated P values after 
correction for multiple testing, controlling the false discovery rate (Holm, 
1 979). Differentially expressed transcripts were detected on the basis of an 
empirical Bayes moderated f-statistic. For the tables shown, values have 
arbitrarily been limited at a nominal two-fold difference in order to present 
the strongest effects only. The Affymetrix CEL, DAT and EXP files have been 
deposited in the Plant Expression database (PLEXdb, accession TA41). SFPs 
were detected between the lines by applying an algorithm described by Xu 
ef al. (2009). 

QTL-specific, constitutive and F. graminearum-respons'we transcripts 
(adjusted P < 0.05) at 8, 24 and 72 hai were calculated by subtracting the 
genetic background provided by NIL4 (susceptible alleles at Fhb1 and 
Qfhs.ifaSA) from NIL2 (Fhb1) and NIL3 (Qfhs-ifa.SA). In addition, we 
generated contrasts from NIL1 and NIL2/NIL3 to yield QTL-specific con- 
trasts in the background of the respective other QTL resistance allele. For 
the identification of transcripts exhibiting constitutive expression changes, 
mock-inoculated samples were compared, whereas, for the F. gramine- 
arum response, transcript abundances of mock samples were subtracted 
from those of F. graminearum-\nocu\ated samples. 

Interactions of the two QTLs were assessed by subtracting abundance 
differences of F. graminearum-responsne transcripts identified in the 
single QTL contrasts (Fhbl, NIL2-NIL4; QfhsJfaSA, NIL3-NIL4) from 



the abundance differences observed in the contrast between the double 
QTL line NIL1 relative to NIL4. Similarly, we calculated the interaction 
effects comparing mock-inoculated samples. The average F. graminearum 
response was calculated as the mean F. graminearum response of all four 
investigated genotypes. 

Annotations, mapping and resequencing of 
Ta.1282.4.S1_at and TaUGT12887 

Annotations were obtained from the wheat HarvEST-database (http:// 
harvest.ucr.edu, Version 1.57). Ontologies, UniProt annotations and EC 
numbers were acquired with the Blast2GO web tool (http:// 
www.blast2go.org). Genes were considered to be 'unknown' if none of the 
search results exhibited E < 10~ 10 . Mapping positions to the identified 
gene transcripts were assigned on the basis of BlastN search results 
against the mapped EST database hosted on the Graingenes website 
(http://wheat.pw.usda.gov/GG2). Moreover, we mapped transcripts 
(>98% identity) to the wheat contig ctg0954b, which harbours Fhbl 
(FN564434.1; Choulet etal., 2010) using GMAP (http://research- 
pub.gene.com/gmap; Wu and Watanabe, 2005). 

Selected candidates were mapped using a collection of Chinese Spring 
nulli-tetra and ditelosomic lines and deletion lines of chromosome 5A 
(Wheat Genetics Resource Center, Kansas State University, Manhattan, KS, 
USA). Ta.1282.4.S1_at and TaUGTI2887 were sequenced from cDNA 
samples, and genomic DNAfrom parental cultivars CM-82036 and Remus. 
All primers are given in Table S10 (see Supporting Information). 

Heterologous expression and functional testing of 
TaUGT12887 in S. cerevisiae 

The coding sequence of TaUGT12887 was optimized for expression in 
5. cerevisiae by reverse translating the amino acid sequence (according 
to 5. cerevisiae codon tables provided by the Kazusa DNA Research Insti- 
tute, Chiba, Japan). The gene was custom synthesized (Eurofins MWG 
Operon, Ebersberg, Germany) including flanking 5' Hind\\\ and 3' Wofl 
restriction sites, which were used to clone the ORF into the Hind////Not/ 
cloning sites of yeast expression vector pYAK7 in frame behind 
the N-terminal c-Myc epitope tag (2\i-P AD m-c-Myc-MCS-T ADm , LEU2; 
Poppenberger etal., 2003). Yeast strain YZGA515 (MATa Ieu2-A7 trpl- 
A63 ura3-52 his3-&200 Iys2-80l ade2-IOi pdr5::TRPI pdr10::hisG 
pdr15:loxP-KANR-loxP aytl::URA) was transformed with the resulting 
construct, pWS60, and also with pWS1921, harbouring the previously 
reported DON-UGT HvUGTI3248 (Schweiger ef al., 201 0) and the empty 
vector (pBP910, Poppenberger etal., 2003), as positive and negative con- 
trols, respectively. 

Crystallized DON was dissolved in water to a stock concentration of 
1 0 g/L. Yeast transformants were selected on synthetic complete medium 
lacking leucine, and grown in liquid selective medium to the exponential 
growth phase [optical density at 600 nm (OD 60 o) = 0.7]. Three microlitres 
of liquid culture dilutions (OD 60 o = 0.05 and 0.005) were spotted on plates 
with yeast extract-peptone-dextrose (YPD) medium containing increasing 
amounts of DON, and incubated for 3 days at 30 °C. Immunodetection was 
performed with a primary mouse anti-c-Myc antibody (1:5000, clone 
9E10; Invitrogen, Carlsbad, CA, USA) and a horseradish peroxidase (HRP)- 
conjugated secondary anti-mouse antibody (Jackson, West Grove, PA, 
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USA). For signal detection, we used the SuperSignal West Pico chemilumi- 
nescent substrate (Pierce, Rockford, IL, USA). 

RT-qPCR validation 

A compiled list of all primers and corresponding transcripts is given in 
Table S1 1 (see Supporting Information). One microgram of total RNA was 
subjected to reverse transcription (QuantiTect Reverse Transcription Kit; 
Qiagen, Hilden, Germany). RT-qPCRs were performed for 45 cycles in 
triplicate using 1 ul of cDNA and primers (300 nM in assay) in combina- 
tion with the MESA GREEN qPCR Mastermix Plus (Eurogentec, Seraing, 
Belgium). The reference genes glyceraldehyde-3-phosphate dehydroge- 
nase (GAPDH), transcription enhancer factor-i (TEF1) and ubiquitin were 
used for normalization. Calculations of the relative fold changes were 
conducted as described in Pfaffl (2001). Transcripts with Ct values >40 
cycles were considered as not present; the corresponding expression ratios 
were displayed as >1 00-fold expressed. 

Metabolite extraction and high-performance liquid 
chromatography-tandem mass spectrometry 
(HPLC-MS/MS) analysis 

Five-hundred milligrams of ground plant material were extracted by the 
addition of 4 mL of solvent consisting of 79% acetonitrile, 1 % acetic acid 
and 20% water. Samples were vigorously shaken at room temperature for 
2 h, spun down and the clear supernatant was subjected to HPLC-MS/MS 
analysis. Fusarium graminearum-treated samples of NILs at 72 hai were 
analysed for DON and D3G (Berthiller etal., 2005). 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online 
version of this article at the publisher's web-site: 

Fig. S1 Reverse transcription-quantitative polymerase chain 
reaction (RT-qPCR) confirmation of selected transcripts. 
Fig. S2 Multiple alignment of the resequenced wheat lipid trans- 
fer proteins (LTPs) corresponding to Ta.1282.4.S1_at in cv. 
CM-82036 and Remus. 

Fig. S3 Immunoblotting of the N-terminal c-Myc antibody 
used to detect uridine diphosphate-glycosyltransferase (UGT) 



MOLECULAR PLANT PATHOLOGY (201 3) 14(8), 772-784 
2013 THE AUTHORS. MOLECULAR PLANT PATHOLOGY PUBLISHED BY JOHN WILEY & SONS LTD AND BSPP 



RNA profiling of Fusarium resistance QTLs 785 



expression levels (A) and the corresponding Coomassie-stained 
acryl amide gel as loading control (B). 
Table S1 Number of transcripts exhibiting a significant difference 
in response to Fusarium graminearum for either Fhb1 or 
Qfhs.ifaSA at the given time point. 

Table S2 Functional classification of Fusarium graminearum- 
responsive transcripts and transcripts showing constitutive quan- 
titative trait locus (QTL)-specific differences in transcript 
abundance (adjusted P < 0.05, >2 fold change). 
Table S3 Transcripts exhibiting a negative difference in response 
for the resistant allele of Qfhs.ifaSA (adjusted P < 0.05, >2 fold 
change, FC) after inoculation with Fusarium graminearum spores 
at 72 h after inoculation. 

Table S4 Transcripts exhibiting constitutive differences in tran- 
script abundance for Qfhs.ifaSA (adjusted P < 0.05, >2 fold 
change, FC). 

Table S5 Transcripts exhibiting constitutive differences in tran- 
script abundance for Fhb1 (adjusted P < 0.05, >2 fold change, 
FC). 

Table S6 Fusarium graminearum-responswe transcripts exhibit- 
ing nonadditive effects between Fhb1 and Qfhs.ifaSA (adjusted 
P < 0.05, >2 fold change, FC). 



Table S7 Transcripts exhibiting nonadditive effects between 
Fhb1 and Qfhs.ifaSA (adjusted P < 0.05, >2 fold change, FC) in 
mock-inoculated samples. 

Table S8 Transcripts exhibiting a significant difference in 
response between NIL1 and NIL3 yielding Fhb1 -related transcripts 
in the presence of the Qfhs.ifaSA resistance allele (adjusted 
P < 0.05, >2 fold change, FC) after inoculation with Fusarium 
graminearum spores. 

Table S9 Transcripts exhibiting a significant difference in 
response between NIL1 and NIL2 yielding Qfhs.ifaSA-related 
transcripts in the presence of the Fhb1 resistance allele (adjusted 
P < 0.05, >2 fold change, FC) after inoculation with Fusarium 
graminearum spores. 

Table S10 Sequences of primers used for sequencing and 
mapping of Ta.1282.4.S1_at and TaUGT12887. 
Table S1 1 Sequences of primers used for the validation of tran- 
script abundances by reverse transcription-quantitative polymer- 
ase chain reaction (RT-qPCR). 
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